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ABSTRACT 

A  model  of  the  radar  evaporative  duct  is  described,  based  on 
Monin-Obukov  similarity  theory.  The  similarity  functions, 
Monin-Obukov  stability  length  and  scaling  parameters  are 
derived  from  experimental  data.  Graphs  are  presented  for 
duct  heights  and  the  number  of  trapped  modes  at  10  GHz  for  a 

wide  range  of  environmental  conditions.  _ _ n 
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1.  INTRODUCTION 

A  decrease  in  refractive  index  with  height  causes  rays  of  electro¬ 
magnetic  radiation  to  bend  towards  the  earth's  surface.  The  mean  rate  of 
change  of  refractive  index  near  the  earth's  surface  determines  the  earth's 
'effective  radius'  and  radar  horizon  [1],  When  the  rate  of  change  of 
refractive  index  equals  a  critical  value,  the  earth's  effective  radius  and 
radar  horizon  become  infinite.  This  is  the  condition  for  super-refraction. 
Greater  rates  of  change  result  in  radar  'ducts',  for  which  the  field 
strengths  for  microwave  propagation  are  not  calculable  by  simple  ray 
theory  and  require  calculation  by  mode  theory  [2]  or  hybrid  models. 

The  input  to  any  such  theoretical  model  is  information  about  s 

(1)  the  vertical  rate  of  change  of  refractive  index  in  the  atmosphere 

and  (in  the  case  of  surface-based  ducts) 

j 

(ii)  the  reflectivity  of  the  rough  sea  surface.  - 

Only  the  latter  is  adequately  described  by  both  theory  and  experimental 
data  [1].  Surface-based  ducts  near  land  may  be  advective  [3],  but  well 
removed  from  continental  effects  and  at  low  altitudes  the  most  important 
mechanism  is  surface  evaporation.  The  refractive  index  variation  for 
elevated  ducts  and  advective  ducts  can  be  obtained  from  radiosonde 
measurements  of  temperature,  humidity  and  height  or  pressure.  In  the  case 
of  surface  evaporative  ducts,  which  are  only  5-25  metres  in  height,  these 
variables  cannot  easily  be  measured  with  sufficient  precision  and  a  profile 
model  must  be  used.  Typically  these  models  require  as  inputs  only  a  few 
simple  bulk  environmental  parameters  such  as  air  and  sea  temperature, 
relative  humidity  and  wind  speed  at  one  height.  From  these,  profiles  of 
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refractive  index  can  be  deduced  which  are  in  reasonable  agreement  with 
direct  measurements  -  especially  when  the  air  and  sea  temperatures  are 
nearly  equal  [4-9].  This  note  briefly  summarizes  existing  models  and 
experimental  data  for  the  evaporative  duct,  and  presents  relations  and 
algorithms  which  are  valid  over  a  wider  range  of  environmental  conditions. 

x  The  outputs  of  the  model  include 

(i)  refractivity  profiles  and  duct  height,  and 

(ii)  the  frequency-dependence  of  trapping  expressed  as  either  the 

number  of  trapped  modes  at  a  given  frequency,  or  the  cut-off 
frequency.  ■  ' y- — -  -  '/'•  '  s  ■  • "  >  1  \ 

c  - 

2.  THEORY 

2.1  Refractivity 

The  refractive  index  (n)  at  microwave  frequencies  is  dependent 

mainly  on  water  vapour  pressure  (e),  temperature  (T)  and  atmospheric 
pressure  (p).  The  refractivity  (N)  is  described  by  the  Debye  formula  [10], 

N  =  106 (n-1 ) 

77.6  (p  +  4810  e)  ^ 

T  T 


with  p  and  e  in  mbar  and  T  in  Kelvin.  The  rate  of  change  of 
refractive  index  with  height  (z)  is  obtained  by  differentiating  equation  1 


dN 

Iz 


{  aN  i  dp  +  r  3N  \  dT  +  r  aN  ) 
1  Ip  J  Hz  1  TT  1  37  1  a?  1 


de 

dz 


(2) 


Evaluating  (2)  at  sea  level  (p0  =  1013.25,  T0  =  288.15  e0  =  10.13)  for  the 
standard  atmosphere  gives 

dN  =  0.269  d£  -1.263  dT  +  4.495  de  (3) 

dz  dz  dz  dz 

The  change  with  pressure  is  fairly  insensitive  to  climatic  variations,  and 
pressure  at  low  altitudes  is  well  described  by  the  barometric  formula 

p  *  Poe‘2/H°  (4) 

where  H0  =  7300  metres.  The  surface  duct  height  is  defined  as  the 

height  corresponding  to  the  minimum  of  the  modified  refractivity 

M  =  N  +  106  z/a  (b) 

where  a  is  the  earth  radius.  That  is,  the  duct  height  corresponds  to  the 
altitude  where 

dN  =  -0.157  m"1  (6) 

dz 

Combining  this  with  equations  3  to  5  leads  to  the  duct-height  relation 
4.495  (de/dz)  -  1.263  (dT/dz)  =  -0.125  (7) 

The  solution  of  eq  (7)  is  single-valued  for  the  evaporative  duct. 
Once  the  humidity  and  temperature  profiles  are  known  it  is  then  a  straight¬ 
forward  matter,  not  only  to  derive  the  duct  height,  but  also  to  determine 
the  M-profile  which  allows  further  analysis  of  the  frequency-dependence  of 
trapping  in  the  duct  [2]. 
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2.2  Neutral  Profiles 

If  the  boundary  layer  adjacent  to  the  sea  surface  is  neutrally 
buoyant,  which  occurs  when  the  temperature  of  the  air  and  sea  are 
approximately  equal,  the  conditions  are  said  to  be  of  neutral  stability. 
The  boundary  layer  can  be  defined  [5]  as  the  layer  in  which  the  fluxes  of 
sensible  heat  (Hs),  latent  heat  (H^)  and  momentum  (-t)  are  constant  (to 
within,  say,  5%).  These  fluxes  can  be  written  as 


HS  - 

pCp  t  u3 

(8) 

hl  = 

p  L 

(9) 

T  = 

-pUlU3 

(10) 

where  q',  t",u'  and  u'  are  fluctuations  in  moisture,  temperature,  downwind 
component  of  wind  and  vertical  component  of  wind  respectively.  The  overbar 
denotes  an  average  over  sufficiently  long  time,  say  one  hour,  for  the 
averages  to  reach  a  statistically  stable  value.  From  Monin-Obukov 
similarity  theory  [b]  there  exist  a  scaling  temperature  (T*),  scaling 
moisture  (q*),  and  friction  velocity  (u*)  which  facilitate  unique 
descriptions  of  non-dimensional  vertical  profiles  (T/T*),  (q/q*)  and  (U/u*) 
for  given  stability  conditions.  The  fluxes  can  then  be  parameterized  as 


Hl  " 

a  K  p  L  u*  q* 

(ID 

HS  = 

a  K  p  0^  u,  T, 

d?' 

T  = 

-p  U*2 

(13) 
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where  K  is  von 

Karman's  constant,  p  is  the  density  of  air. 

Cp  is  the 

heat  capacity  of 

air,  and  L  is  the  latent 

heat  of 

vapori zation 

.  The 

*  * 

constant  a  is 

the  ratio  of  the  eddy  di f fusi vities 

of  heat  and 

momentum 

which  is  usual ly 

assumed  to  be  unity  [5]  but 

has  been  shown  from 

experiment 

‘•V- 

to  be  around  1.3b  [11].  This  uncertainty  will  not 

affect  the  conclusions 

«  •*.  1 
*  m  m 

in  this  section 

since  we  define  the  scaling 

parameters  according 

to  the 

gradients  under 

neutral  conditions 

■- .  * 

dU  = 

u* 

(14) 

*-•  *  * 

dz 

Kz 

• 

dT  = 

T* 

(lb) 

dz 

Kz 

d^  = 

q* 

(16) 

dz 

Kz 

Integration  of  (14)  -  (16)  leads  to 

the  mean  quantities 

V 

u 

u*  £rt{z/z0) 

(17) 

■  < 

1 T 

T 

T  +  T*  *n(z/z0) 

°  K 

(18) 

*  *  v 

^  | 

q 

qn  +  ll  ^n(z/z0) 
u  < 

(19) 

where  T0  is  the  temperature  and  qQ  is  the  mixing  ratio  at  z=Zq.  The 
latter  is  the  roughness  length  for  momentum  which  has  been  experimentally 
related  to  the  friction  velocity  [12] 


a  u*/g  +  (v/ 


i  a-5.5K 
e 


(20) 


2 

where  a  is  the  Charnock-El 1 ison  constant  («  0.011),  g  =  9.8  m/sec  and 
v  is  the  kinematic  viscosity  of  air  (v  «  14  x  10"6  m2/sec).  The  first 
term  in  eq  (20)  represents  the  increase  in  surface  roughness  with 
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increasing  wind  stress  owing  to  surface  waves,  while  the  second  term 
represents  an  aerodynami cal ly  smooth  surface  and  is  important  for  wind 
speeds  less  than  a  few  metres  per  second.  Typical  values  for  moderate  wind 
(U  =  lOm/sec)  are  u*  =  0.36  m/sec  and  z0  =  1.5  x  10~4  m.  In  some 
models  of  the  evaporative  duct  [6,7]  zq  is  taken  as  a  constant  roughness 
length  with  a  value  of  ~  1.5  x  10*4m.  This  assumption  is  probably 
adequate  as 

(i)  we  are  not  concerned  with  the  behaviour  at  heights  which  are  less 
than  the  wavelength  of  microwave  radiation,  and 

(ii)  the  integration  constants  Tq  and  qQ  are  probably  not 
measurably  different  from  the  bounds  determined  by  the  sea  surface 
temperature. 

The  scaling  factors  for  neutral  profiles  can  then  be  determined 
from  the  expressions 

u*  =  KUj  (21) 

^nTz;/z0) 

T*  *  u*  (T1  -  Tq)  (22) 


-  u*(qj  -  qQ)  (23) 


where  T, ,  u,,  q,  are  the  temperature,  wind  speed  and  humidity  measured  at 
z  =  z;  •  ^or  moderate  winds  u*/y  _  0.U38  m/sec  [4,8]  so  that,  when  T^  Tp 
the  above  relations  reduce,  after  some  substitution,  to  a  simple  expression 
for  the  duct  height  in  terms  of  10-metre  bulk  differences 


E 


with  e  in  mbar  and  T  in  Kelvin.  For  typical  moisture  gradients  of  up 
to  10  mbar  in  10  metres,  neutral  duct  heights  are  up  to  30  metres,  with  the 
temperature  term  being  of  little  importance. 

The  assumption  of  a  constant  drag  coefficient  (u^/U  )  for  a  given 
measurement  height,  and  a  simple  relation  between  wind  and  humidity 
profiles,  leads  to  the  conclusion  that  neutral  duct  heights  are  independent 
of  wind  speed  for  given  humidity  difference.  This  is  a  common  feature  of 
most  of  the  models  in  the  literature  for  neutral  duct  heights  [4,  6-8].  If 
the  magnitude  of  thp  air-sea  temperature  difference  is  significant,  the 
above  models  are  not  adequate  and  the  effect  of  stability  must  be  taken 
into  account. 


Stable  and  Unstable  Profiles 


When  the  air  is  warmer  than  the  sea  the  conditions  are  saiu  to  be 
stable,  and  a  warmer  sea  is  said  to  be  unstable.  According  to  Monin-Obukov 
similarity  theory,  the  gradients  at  height  z  in  the  boundary  layer  can  be 
determined  from  T*t  q*,  u*  and  a  dimensionless  height  (z/L).  The 
calculation  of  the  Monin-Obukov  stability  length  (L)  will  be  outlined  in 
section  2.5. 


The  gradients  of  U,  T  and  q  can  be  written  in  terms  of  stability 


functions  4>(z/L) 
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dU  = 

u* 

♦U(z/L) 

(25) 

clz 

YZ 

d]_  = 

T* 

4>t(z/L) 

(26) 

dz 

Yz 

d£  = 

q* 

♦q(z/L) 

(27) 

dz 

Yz 

Universal  profiles  4^  (z/L)  are  related  to  ( z /L )  [5]  as 

*  (z/L)  =  [  /  1~  *i(z/L)  d(z/L)  (28) 

!  z0/L  z/L 


or 


*1  (z/L)  =  1  -  z  (d*Vdz  ) 


(29) 


This  leads  to  expressions  for  mean  quantities  analogous  to  eqns 
(17)  -  (19), 


U 

=  u* 
K- 

[  Jin 

(z/z0)  -  *u 

(z/L)  ] 

(30) 

T 

=  T 

+  T* 

[  Jin  (z/z0) 

-  *T  (z/L)  ] 

(31) 

0 

K 

q 

=  % 

+  9* 

[  Jin  (z/z0) 

-  Tq  (z/L)  ] 

(32) 

It  follows  from  (28)  that  Tq  and  qQ  are  defined  as  the  values 
of  T  and  q  evaluated  at  z=z0  (the  roughness  length  for  momentum). 
The  error  in  equating  these  to  the  equilibrium  sea-level  values  should  be 
negligible.  A  more  detailed  discussion  is  given  in  Annex  A. 


2.4 


Empirical  Stability  Functions 


All  profiles  have  the  feature  that  as  z/L  +  0,  ^  +  1  and 

'h  +  0.  Under  stable  conditions,  the  stability  functions  can  be 
approximated  by  the  first  two  terms  in  the  Taylor  series  expansion 

$-j  «  1  +  a.j(z/L)  z/L  >  0  (33) 

where  a  =  4.5  [7]  and  ay  =  a  =  6.35  [11]. 

Under  unstable  conditions  the  stability  functions  are  very  non-linear  in 
z/L.  Under  these  conditions  the  stability  function  for  wind  follows  the 
Keyps  profile  [5]  which  can  be  approximated  as  [11] 

4>u  (1  -  18  (z/L)  )~1/4  z/L  <  0  (34) 

Temperature  and  humidity  functions  for  unstable  conditons  are  in  good 
agreement  with  the  expression  [11] 

<i»T  =  <J>q  =  (1-9  (z/L)  )  z/L  <  0  (35) 

From  equation  (28)  the  universal  profiles  follow  as 


The  function  f  (z/L)  for  unstable  conditions  is  calculated  from  a 
fit  of  the  numerical  integration  [5,  7]  of  <j> 
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2.b  Monin  -  Obukov  Stability  Length 


The  dimensionless  quantity  (z/L)  is  the  ratio  of  convective  to 
mechanical  energy  production  in  near  neutral  air  [5],  The  Monin-Obukov 
stability  length  (L )  can  then  be  written  as 

3 

L  =  u*  j  (3‘ 


K  g  (  t'u^  +  0.608  T  q ^ ) 


The  sensible  heat  flux  (p  C t'u'  )  and  moisture  flux  (p  q'u*) 
are  not  easily  measured,  and  various  parameterization  schemes  have  been 
proposed.  Most  are  of  the  form  [13] 


A  +  B  U  aT 


CE  U2  Aq 


where  Uz  is  the  mean  wind  speed  at  the  reference  height,  and  aT  and  Aq 
are  the  air-sea  temperature  and  moisture  differences  respectively.  These 
will  not  be  used  in  this  model  for  reasons  of  internal  consistency,  but  are 
useful  for  parametric  schemes  of  evaporative  duct  height  in  near-neutral 
conditions  [4]. 


A  quantity  which  is  closely  related  to  z/L  is  the  gradient 
Richardson  number  Ri . 


q  (  aT/?z  ) 

T  (  3U/3Z  )2 
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From  the  profile  relations  in  the  previous  sections,  (42)  can  be  re-written 
as 

Ri  =  gz t\l  .  f  (in  (z\/z0)  -  yu)2  .  6T  (z/L)  •,  (43) 

T  U2  (*n  (z  1  /z0 )  -  H-t)  ^  (z/L) 

For  near-neutral  conditions,  with  L)  and  T  measured  at  6  metres,  the  term 
in  brackets  is  of  order  10.  Under  these  conditions,  (43)  reduces  to  the 
definition  of  the  bulk  Richardson  number  used  in  refs  [6]  and  [7] 

R  =  IP  .3“!.  (44) 

,b  T  U2 

In  ref. [7],  R . b  is  first  calculated  according  to  eqn  (44)  and  z/L  is  then 
calculated  by  a  simple  empirical  fit. 


In  this  model  we  prefer  to  use  eqn  39  directly,  with  terms  in 
t'u'  and  q'u'  replaced  by  terms  in  scaling  temperatures  and  moisture 
respectively 


L 


2 

u*  T 


1.35  K  g  (T* 


0.608  T  q*) 


(45) 


where  the  factor  1.35  arises  because  T*  and  q*  are  defined  here  in  terms 
of  the  respective  gradients  and  not  the  fluxes.  Businger  [11]  has  measured 
this  quantity  which  is  a  ratio  of  eddy  diffusivities,  by  comparing  the  two 
methods  of  calculation  from  flux  and  profile  measurements.  The  difficulty 
with  (45)  is  that  U*,  T*  and  q*  are  themselves  functions  of  L 


u*  =  KUj/  r*n  (Zj/z0)  -  (Zj/L)  1 


(46) 


K  (T  j-  TQ)  /  ' in  (z  x/z  Q )  -  tt  (Zj/L)  ' 


K  (qt-  qQ)  /  ( in  (z  j/zo )  -  v  (Zj/L)  1 


The  model  solves  (45)  -  (48),  which  are  well-behaved  functions,  by 
a  stepwise  refinement  iterative  procedure.  Given  a  reasonable  first  quess 
for  L,  the  test  quantity  1/L  converges  to  10"4  metre  typically  in  3-4 
iterations. 


Comparisons  of  R-^,  z^/i  (ref. 7)  and  z 1 /L  (iterative  solution) 
are  shown  in  figure  1.  The  bulk  Richardson  number  as  defined  in  (42)  is 
not  a  particularly  good  estimator,  since  the  density  gradients  are 
approximated  by  temperature  gradients,  and  does  not  account  for  the  fact 
that  moist  air  is  less  dense  than  dry  air.  According  to  similarity  theory 
the  Richardson  number  is  a  unique  function  of  L  and  so  figure  1  could  be 
used  as  a  look-up  graph.  However,  in  order  to  compare  L  from  either  a 
parametric  model  (eqns  39-41)  or  equation  45,  the  Richardson  numbers  need 
to  be  defined  in  terms  of  virtual  temperature  (T  ) 


T  (  1  +  0.608  q  ) 


aT  +  0.107  Ae 


(49(a)) 


in  which  q  is  the  mixing  ratio  and  Ae  is  the  change  in  water  vapour 
pressure  in  mbar.  Typical  Ae  for  tropical  waters  is  -10  mbar  [4],  so 
that  neutral  conditions  occur  at  aT  »  +1  deg.  Celsius. 


From  figure  1  it  is  clear  that,  over  a  moderate  range  of 


temperatures  (and  therefore  saturated  vapour  pressures)  and  relative 
humidities,  zx/l  calculated  by  the  method  of  ref.  7  is  not  sufficiently 
accurate  for  the  model  described  here.  It  is  used,  however,  as  a  first 
guess  for  the  iterative  algorithm  for  L,  except  for  weakly  unstable 
conditions  where  it  predicts  the  wrong  sign  for  L. 


2.6  Calculation  of  the  Evaporative  Duct  Height 

The  evaporative  duct  height  is  defined  as  the  height  at  which  the 
modified  refractivity  (M)  is  a  minimum.  From  eq  (7), 

-0.125  =  4.495  c  ,  d£  >  -1.263  (  dT  )  (50) 

dz  ‘  dz 

where  c  is  the  ratio  of  water  vapour  pressure  in  mbar  to  mixing  ratio 
(  =  1630).  It  is  convenient  here  to  define  a  new  variable  N',  closely 
related  to  the  refractivity,  for  which  the  rate  of  change  of  N'  with 
height  is  the  rate  of  change  of  N  less  that  due  to  pressure  alone; 

dff  =  4.495  cf  d£  1  -  1.263  ,  dT  n  (51) 

dz  1  dz  1  1  dz  j 

4.495  c  q*  4>q(z/L )  -  ^63  T«  *  (Z/L)  (52) 

Kz  Kz 

Now  if  $  =  4>q  ,  N'  follows  the  similarity  rules 


dN' 

dz 

where  4>N 
K 


N*  (z/l) 

Kz 

=  <j>q  and 

4.495  c  q,  1.263  T* 
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The  duct-height  (z)  is  then  obtained  from  the  relation 

N*  6n(z/L)  =  -0.12b  K  z  (5b) 

In  the  program,  equation  (55)  is  solved  by  the  secant  iterative  method,  and 
z  usually  converges  to  with  10~3  metre  after  2  iterations. 

Results  for  two  extremes  of  relative  humidity  (30%  and  90%)  are 
shown  in  figures  2-6  for  Ts  =  15,  18,  21,  24  and  28  degrees  Celsius.  In 
all  cases  duct  height  increases  with  wind  speed  in  the  unstable  region 
(Ta  -  Ts  <  0)  and  is  wind-speed-independent  for  neutral  conditions.  The 
effect  is  reversed  for  stable  conditions. 

Results  for  the  lowest  wind  speed  (3  m/sec)  are  not  plotted  past 
aT  =  0.5°C  since  this  corresponds  to  z/L  >  0.1,  which  is  the  limit  of 
validity  of  this  model  and  of  similar  models  based  on  turbulence  theory  [5]. 
The  lower  limit  (z/L  =  -1)  seems  a  less  important  bound  since  the  scaling 
laws  extrapolate  adequately  in  the  unstable  region  [8], 


2.7  Comparison  with  other  Models 

Most  models  of  the  surface  evaporative  duct,  which  can  be  found  in 
the  literature,  provide  similar  results  for  neutral  conditions.  This  is 
because  the  neutral  drag  coefficient  is  independent  of  wind  speed,  so  that 
the  duct  height  is  approximately  linear  in  the  humidity  difference.  The 
results  of  Jones  and  Stewart  [4]  reduce,  after  some  substitution  to 

d  =  2.50  (eQ-  ejn)  (56) 


lb 


where  e1Q  is  the  water  vapour  pressure  in  mbar  at  10  metres,  and  eQ  is  the 
sea  level  value.  The  difference  in  the  numerical  factor  between  this 
result  and  equation  (24)  is  due  to  different  values  used  for  the  friction 
velocity  and  heat  transfer  coefficients.  The  Hitney  model  [7]  (also  used 
in  the  NOSC  ‘  I  RE  PS '  computer  program)  covers  unstable,  neutral  and  stable 
conditions  and  reduces  to  a  closed  form  expression 


d 


_ an _ 

bj  (  An(z1/z0)  +  5.2  Zj/L  )  -  5.2  aN/L 


(57) 


for  stable  conditions  and  an  analogous  expression  for  unstable  conditions. 
Although  no  derivation  is  given,  the  model  is  similar  to  that  used  by 
Jeske  [6]  which  assumes  identical  stability  functions  for  wind,  temperature 
and  moisture.  It  also  appears  to  assume  that  9N/3z  is  linear  combination 
of  Ae  and  AT.  Under  neutral  conditions  the  Hitney/Jeske  models  reduce  to 


d  =  0.025  +  3.4  (eQ  -  e1Q)  (58) 

which  is  very  close  to  the  result  of  the  model  (eq  24).  For  stable  and 
unstable  conditions,  the  general  trends  with  temperature  difference  and 
wind  speed  are  similar,  the  main  difference  being  due  to  the  use  of 
normalized  wind  profiles  for  temperature  and  moisture. 

The  USN  Weather  Service  model  [8]  is  closer  to  that  described  here 
except  that,  in  order  to  avoid  iterative  solutions,  the  equality  of  wind, 
temperature  and  moisture  stability  functions  was  also  assumed.  The  USN 
model  reduces  to  the  genera]  expressions 

d  -  *u  (Zl/L)  (A  q*  -  B  T*  ) 


where  A  and  B  are  constants  determined  by  the  selected  reference 
1  evel  (z  ). 


(59) 


2.8 


Frequency-dependence  of  Duct  Strength 


The  ability  of  a  duct  to  support  trapped  modes  increases  with  duct 
height,  frequency  and  aM  [2].  In  order  to  calculate  the  signal  strength  in 
the  duct  it  is  necessary  to  use  mode  theory  [2,9],  but  this  is  not  feasible 
for  routine  environmental  predictions. 

The  phase  criteria  for  trapped  rays  can  be  readily  calculated  from 
ray  theory  by  considering  the  ray  to  be  composed  of  two  components  -  one 
ray  propagating  along  the  duct  axis,  and  a  standing  wave  in  the  vertical 
direction.  The  condition  for  maximum  energy  flow  along  the  duct  requires 

that  ‘similar  points'  for  a  ray  be  in  phase,  and  this  is  calculable  by 

considering  only  the  vertical  ray  and  two  reflections  -  one  at  the  sea 
surface  and  the  other  at  some  turning  point,  bm  [2]>  The  graz1ng  angle  for 

the  surface  reflection  is  sufficiently  small  for  the  phase  change  on 

reflection  to  be  equated  to  180°,  and  so  the  phase  condition  is 

2  rr x  10^  I  m  /  M ( z ) -M ( d )  dz  =  m  "  1/4,  m  =  1,2,3..  (60) 

\  o 

where  d  is  the  duct  height.  When  bm  is  set  to  d,  the  solution  of 
(60)  yields  the  number  of  trapped  modes  (m).  When  m=l,  the  solution  gives 
the  maximum  trapped  wavelength  \  or  minimum  trapped  frequency,  or  cut-off 
frequency,  (fCQ). 

In  general,  equation  (60)  cannot  be  solved  analytically  and  the 
M-profile  must  be  broken  up  into  integrable  segments.  In  the  model  n  =  100 
linear  segments  are  used,  but  n  as  low  as  10  may  give  adequate  results 
for  most  purposes  as  long  as  the  M-profile  is  well  described  in  the  lowest 
10  or  so  metres,  where  most  of  the  M-deficit  occurs. 


In  each  linear  segment  a  slope  ( k ^ )  and  intercept  (M  ^  )  are 
defined  so  that,  in  the  ith  segment, 

*yz)  =  M0i  -  kiz  (6D 

and  equation  (6U)  can  be  rewritten  as 

_  n-1  (i+l)d/n 

m  -1/4  =  2  A?  x  10  l  j 

x  i=o  id/n 

■  2  CL »  »°~3  I  -  2  k,1*  d3^  [  (1  -  m)3'2  -  (1  -1J3/2  ]  (62) 

X  i  =0  3  n  n 

where 

ki  =  "(Mi+l  -  Mi)  /  d  (63) 

for  equal  sized  segments. 

Setting  m  =  1  and  solving  for  frequency  (c/x)  yields  the  nominal 
cut-off  frequency  (fCQ),  which  is  one  of  the  outputs  of  the  program.  The 
number  of  trapped  modes  at  any  other  frequency  (f)  is  then  simply 

m  =  0.2b  +  0.7b  (f/f  )  (64) 

Equation  (64)  allows  non-integer  values  of  m.  In  ray  theory  the  turning 
point  is  well-defined  for  any  ray  so  that  only  the  integer  part  of  m  is 
significant.  (That  is,  m  <  1  implies  no  trapping).  The  interpretation 
^rom  mode  theory  is  that  m  is  in  the  order  of  the  number  of  trapped  modes 
that  need  to  be  calculated  for  the  contribution  from  trapped  modes  to 
approach  convergence. 

The  calculated  modified  refractivity  profiles,  M(z),  for  several 
air-sea  temperature  differences  are  plotted  in  figure  7.  The  dashed  line 


is  the  "typical"  refractivity  profile  cited  in  reference  [2]  which 
facilitates  a  rough  estimate  of  cut-off  frequency  based  only  on  the  duct 
height  and  m  =  1.  Inspection  of  equation  (60)  and  figure  7  suggests  that 
this  typical  profile  may  over-estimate  or  under-estimate  the  cut-off 
frequency  depending  on  stability.  In  the  case  of  super-refraction  (M 
fairly  constant)  the  results  would  be  spurious. 

The  general  shape  of  the  profiles  in  figure  7  also  suggests  that: 

a.  experimental  measurement  of  duct  height  would  normally  be  highly 
inaccurate  since,  by  definition,  zM/?z  =  0  at  the  duct  height, 
and 

b.  with  equation  (60)  it  indicates  that  significant  variations  in  the 
limit  of  integration  (the  estimated  duct  height)  nay  have  little 
effect  on  the  computed  value  of  m,  as  long  as  M  is  well  described 
in  the  lowest  10  or  so  metres. 

The  effect  of  stability  on  the  number  of  trapped  modes  at  10  GHz 
is  shown  in  figures  8-12.  The  trends  are  similar  to  those  of  the  duct 
height  -  namely  that  the  number  of  trapped  modes  increases  with  wind 
strength  when  the  sea  is  warmer  than  the  air  and  vice  versa.  Figs  8-12 
also  indicate  that  when  the  humidity  is  high,  only  one  or  two  modes  will  be 
trapped  at  10  GHz.  This  conclusion,  as  well  as  the  dependence  of  duct 
strength  on  stability,  was  observed  in  very  early  measurements  of  the 
evaporative  duct  [14,15]. 


2.9 


Propagation  Modelling 


I 

i 

t  Other  outputs  of  the  program  include  the  M-profile  between  sea 

! 

level  and  100  metres,  in  1-metre  segments.  This  output,  together  with 
standard  models  for  the  rest  of  the  tropoJenere  [1],  is  the  required  input 
into  microwave  propagation  models  whether  based  on  mode  theory,  ray  theory 
or  a  hybrid  approach  [2,9]. 

The  simplest  description  of  the  distribution  of  electromagnetic 
energy  with  range  and  height  is  the  ray  trace.  Figure  13  was  obtained  by 
piping  the  output  of  the  evaporative  duct  program  into  a  simple  ray  tracing 
program.  Work  is  in  progress  to  incorporate  the  evaporative  duct  model 
into  multipath  propagation  programs  [1,17]  to  produce  vertical  coverage 
diagrams  under  ducting  conditions. 


3. 


SUMMARY 


3.1  The  model  described  in  this  note  calculates  duct  height  and  an 
estimate  of  the  number  of  trapped  modes,  for  the  surface  evaporative  duct 
under  unstable,  neutral  and  slightly  stable  conditions. 

3.2  The  empirical  inputs  into  the  model  are  temperature  and  humidity 
stability  functions.  Refractivity  profiles  are  uniquely  determined  by 
these  inputs  but  have  not  been  directly  measured  with  sufficient  precision 
over  a  range  of  stabilities  to  compare  with  the  output  of  the  model. 

3.3  No  account  is  made  of  the  behaviour  at  high  wind  speeds,  where  the 
effects  of  spray  and  breaking  wave  phenomena  [16]  may  need  to  be 
incorporated. 

3.4  In  the  case  of  neutral  stability,  the  model  gives  similar  results 
to  other  models  in  the  literature.  For  stable  and  unstable  conditions,  no 
simple  analytical  expressions  exist  for  the  duct  parameters,  and  iterative 
or  recursive  algorithms  must  be  used.  Representative  calculations  show 
that  the  convergence  of  the  algorithms  used  in  the  model  is  very  rapid. 

3.5  The  output  of  the  associated  computer  program  includes  the 
modified  refractivity  profile,  which  is  the  required  input  for  microwave 
propagation  models. 
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ANNEX  A  TO 


RANRL  TECH  NOTE  3/85 
DTTTTD  APR  85 


ROUGHNESS  LENGTH  FOR  MOMENTUM,  WATER  VAPOUR  AND  SENSIBLE  HEAT 


If  measurements  of  wind,  humidity  and  temperature  are  made  at  two 
heights,  the  friction  velocity  (u*),  scaling  moisture  (q* )  and  scaling 
temperature  (T*)  can  be  obtained  from  the  expressions 


U* 

K  (U2  -  Uj) 

(A-l) 

*n(z,/Zj)  -  VS71) 

+  WL) 

q* 

k  (q 2  ‘ 

( A -2 ) 

en (z  ?/z  i )  -  (z  p/L  ^ 

+  "q(z/L) 

t*  = 

K  (T2  -  Tj) 

( A  -  3 ) 

Jn  (Zj/Zj)  -  ,1'T  ( 2  ?  /L  )  +  v  (Zj/L) 


where  the  functions  ^(z/L)  are  obtained  empirically  as  described  in 
section  2.4.  Routine  measurements  at  sea  are  made  at  only  one  height, 
usually  around  10  metres  above  sea  level.  The  second  reference  level  is 
often  the  roughness  length,  or  an  estimate  of  it.  The  roughness  length  for 
momentum  (z0)  is  the  height  at  which  the  wind  extrapolates  (using  eqns  30 
and  36)  to  zero.  Similarly  the  roughness  length  for  sensible  heat  ( z-j- )  is 
the  height  at  which  the  temperature  extrapolates  to  the  sea  surface 
temperature  (T  ),  and  the  roughness  length  for  water  vapour  (zq)  is  the 
height  at  which  the  absolute  humidity  corresponds  to  the  saturated  humidity 
above  sea  water  (q  ).  It  is  a  common  feature  of  several  models  of  the 
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evaporative  duct  to  assume  that  Zy  =  zq  =  z0,  and  it  is  also  assumed  in 
some  models  that  z0  has  a  constant  value  of  1.5  x  10‘4  m  [6,7].  From 
eqn  (20),  z0  increases  linearly  with  wind  speed  for  moderate  wind  speeds, 
so  that  the  latter  assumption  is  not  valid.  The  roughness  lengths  for 
water  vapour  and  sensible  heat  may  also  differ  significantly  from  z0  Since 
momemtum  is  transferred  across  the  air-sea  interface  by  both  molecular 
diffusion  and  pressure  forces,  while  heat  and  mass  transfer  are  essentially 
controlled  by  molecular  diffusion  alone.  For  moderate  wind  speeds  z0  is 
generally  greater  than  Zy  and  zq  since,  in  the  case  of  momentum,  the 
roughness  elements  generate  drag  due  to  pressure  forces,  whereas  in  the 
case  of  admixture,  these  same  elements  cause  a  sheltering  effect, 
inhibiting  transfer  at  the  surface  [18]. 

The  ratio  of  Zy  and  z^  to  zQ  can  be  calculated  from  the 
expressions  [18] 

V zQ  =  7.4  exp(  -7.3  K  Sc 1/2  R.l/4  )  (A-4) 

Zy/Z0  =  7.4  exp(  -7.3  K  Pr1/2  )  (A-5) 

where  Sc  is  the  Schmidt  number,  or  ratio  of  kinematic  viscosity  of  air 
(v)  to  molecular  diffusivity  of  water  in  air  (»  0.595), 

Pr  is  the  analogous  Prandtl  number  (>»  0.71)  for  heat  diffusivity, 

K  is  von  Karman's  constant  {»  0.4),  and  RQ  is  the  roughness  Reynolds 

number 

Ro  =  zo  u*/v  (A_6) 

Comparing  eqs  (A-6)  and  (20),  R  reduces  for  moderate  wind  speed  to 
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R  =  8U. 2  u 

o  * 


with  u+  in  m/sec  so  that 


(A-7) 


2j/z0  =  7.4  exp{-7.4  u*  ) 


zq/2o 


7.4  exp(-6.  7  ujf2) 


(A-8) 


(A-9) 


For  a  moderate  wind  speed  of  Uio  =  lOm/sec  (u*  =  0.38  m/sec), 

zj/zQ  =  0.079  and  zq/z0  =  0.116.  For  near-neutral  conditions,  the  error 

in  evaluating  the  scaling  temperature  and  moisture  (T*,  q*)  from  eqns  (18) 
and  (19)  is  then  approximately  20%  in  both  cases. 

Under  near-neutral  conditions,  the  evaluation  of  T  and  q  at  some 
height  (z?)  other  than  the  measurement  height  (z^)  can  be  obtained  from 

eqns  (A-2)  and  (A-3)  with  ?(z/L)  =  0,  without  explicity  calculating  T*  and 


T0  +  "  l”o )  tn  (z 


q0  +  (qi  -  qo)  *n  (VV 


(A-10) 


(A-ll) 


An  order-of-magnitude  error  in  Zy  and  z^  would  therefore  propagate 
an  error  of  only  a  few  percent  of  the  corresponding  air-sea  difference  for 
heights  of  1  to  100  metres  when  the  reference  height  is  10  metres. 
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Similar  conclusions  apply  for  non-neutral  stability  where 
differences  between  model  predictions  are  more  likely  to  be  dominated  by 
differences  in  stability  functions  <j>(z/L)  and  ^(z/L).  It  is  not  clear, 
however,  that  the  stability  functions  derived  from  flux  profile 
measurements  are  themselves  not  biased  by  the  choice  of  Zy  [11]. 

There  is  generally  a  paucity  of  experimental  data  for  simultaneous 
moisture,  heat  and  momentum  flux  measurements  [19],  which  are  required  for 
the  separation  of  the  variables  described  in  this  Annex.  Estimates  of 
similarity  functions,  roughness  lengths,  scaling  temperature  and  moisture, 
and  von  Karman's  constant  from  limited  or  independent  experiments  are  all 
susceptible  to  the  type  of  error  described  above. 
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Figure  3.  Calculated  evaporative  duct  height  for  sea  temperature  of  18  deg  C. 
Wind  speed  *  3,  6,  9,  12,  and  15  m/sec 
Relative  humidity  t  302  and  902 
Measurement  height  i  10  metres 
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Figure  4.  Calculated  evaporative  duct  height  for  sea  temperature  of  21  deg  C. 
Kind  speed  :  3,  6,  9,  12,  and  15  m/sec 
Relative  humidity  i  3CX  and  90Z 
Measurement  height  t  10  metres 


Air-sea  temperature  difference  (  deg  Celsius  ) 

Figure  5.  Calculated  evaporative  duct  height  for  sea  temperature  of  24  deg  C. 
Wind  speed  i  3,  6,  9,  12,  and  15  m/sec 
Relative  humidity  i  30X  and  90Z 
Measurement  height  i  10  metres 


Air-sea  temperature  difference  (  deg  Celsius  ) 

Figure  6.  Calculated  evaporative  duct  height  for  sea  temperature  of  28  deg  C 
Kind  speed  *  3,  6,  9,  12,  and  15  m/'sec 
Relative  humidity  i  302  and  902 
Measurement  height  t  10  metres 


Height  above  Sea  Level  (  metres  ) 


Modified  Refractivity  (  M  units  ) 

Figure  7.  Modified  refractivity  for  7  air-sea  temperature  differences. 
Sea  temperature  i  28  deg  Celsius 
Kind  speed  t  6  m/sec 
Relative  humidity  j  85Z 
Measurement  height  i  10  metres 
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Numbar  of  Ircppad  aodaa  OO  at  10  GHz  for  wo  ta^w-atura  of  15  dag  C. 
Wind  speed  i  3a  6,  9a  12a  and  15  «/sec 
Relative  humidity  i  XX  and  9QX 
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